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ABSTRACT 

The colossal power output of active galactic nuclei (AGN) is believed to be 
fueled by the accretion of matter onto a supermassive black hole. This central 
accreting region of AGN has hitherto been spatially unresolved and its structure 
therefore unknown. Here we propose that a previously reported 'deep minimum' 
in the X-ray intensity of the AGN MCG— 6— 30— 15, was due to a unique X-ray 
occultation event and that it probes structure of the central engine on scales 
< 10 14 cm, or 1.4 x 1CT 7 arcseconds. This resolution is more than a factor of 
~ 3 x 10 6 greater than is possible with current X-ray optics. The data are 
consistent with a bright central source surrounded by a less intense ring, which 
we identify with the inner edge of an accretion disk. These may be the first direct 
measurements of the spatial structure and geometry of the accreting black-hole 
system in an active galaxy. We estimate a mass lower limit for sub-Eddington 
accretion of 3.1 x 1O 5 M . If the ring of X-ray emission is identified with the 
inner edge of an accretion disk, we get mass upper limits of 1.9 x 10 8 and 
9.1 x 10 8 M Q for a non-rotating and maximally rotating black hole respectively. 
We point out that our occultation interpretation is controversial in the sense 
that X-ray variability in AGNs is normally attributed to intrinsic physical 
changes in the X-ray emission region, such as disk or coronal instabilities. 

Subject headings: accretion disks - black hole physics - galaxies:active - 
galaxies: individual: MCG —6—30—15 - X-rays: galaxies 
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1. Introduction 

The accretion of matter onto a supermassive black hole as a mechanism for fueling 
the output of active galactic nuclei (AGN) is a paradigm strongly supported by recent 
spectroscopic observations of the iron Ka X-ray emission line (Tanaka et al. 1995; Yaqoob 
et al. 1995; Nandra et al. 1997a and references therein). The extreme Doppler and 
gravitational energy shifts of the line photons, together with the shape of the line are 
consistent with an origin in a disk rotating about a black hole (Fabian et al. 1989; 
Laor 1991). Strong gravitational redshifts, in which photon energies are changed by 
more than 10%, occur only when matter approaches closer than ~ 20 gravitational radii 
(= 20r g ;r g = GM/c 2 ) from a compact object. However, it is not possible to directly 
map the physical structure of the system since the highest spatial resolution of X-ray 
optics technology is a factor ~ 10 6 too poor for even the closest AGN. Optical and radio 
observations provide greater resolution but the bulk of the emission at these wavelengths 
is not generated close enough to the central engine. So far, the highest spatial resolution 
observations, at radio wavelengths, have revealed a Keplerian disk in the AGN NGC 4258 
down to only ~ 60, 000r 9 (Miyoshi et al. 1995; Maoz 1995). This still falls short by a factor 
~ 3000 of mapping the black-hole region. 

The AGN MCG-6-30-15 (z = 0.008) was observed by the X-ray astronomy satellite 
ASCA {Advanced Satellite for Astrophysics and Cosmology; Tanaka, Inoue, and Holt 1994) 
for rsj 4.2 days on 23 July 1994. Results from this observation have already appeared in the 
literature, including the 0.5-10 keV lightcurve (Iwasawa et al. 1996, hereafter 196; Reynolds 
1997; Yaqoob et al. 1997; see Figure 1) and a broad, asymmetric, variable iron K line with 
a strong red wing, consistent with a disk inclined at ~ 30° rotating about a black hole 
(Tanaka et al. 1995; 196). The X-ray luminosity exhibits erratic variability on all timescales 
down to < 50 s (Matsuoka et al. 1990; Green, McHardy, and Lehto 1995; Reynolds et al. 
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1995; Nandra et al. 1997b). Causality arguments alone cannot put constraints on the size 
of the X-ray emission region since the high-frequency, lower amplitude variability may occur 
at localized regions of the source. Figure 1 shows an extended intensity dip at the end of 
the observation, from ~ 3.3 x 10 5 s to ~ 3.6 x 10 5 s. This feature has previously been dubbed 
as the 'deep minimum', or DM (196). A closer inspection (Figure 2a) reveals a remarkable 
(albeit approximate) symmetry about the minimum luminosity. 

We propose that the dip was caused by an occultation of the X-ray source by 
optically-thick matter. This interpretation is controversial and 'non-standard', as X-ray 
variability in AGNs is normally attributed to intrinsic properties of the X-ray emission 
region, such as disk or coronal instabilities. However, so little is know about structure of 
the central engine in AGNs that the occultation scenario should be explored. The obscurer 
must be optically thick because the dip continuum spectrum only shows evidence for weak 
absorption, nowhere near enough to explain the observed intensity variation over the whole 
ASCA bandpass (see Weaver & Yaqoob 1998, hereafter WY98). The luminosity at the 
absolute minimum of the dip is ~ 0.4 of the pre-dip value and must represent persistent 
emission which has much smaller surface brightness than the primary source. Hereafter we 
will refer only to the primary X-ray emission, unless explicitly referring to the persistent 
emission. The proposed obscurer very likely hides the most compact and variable part 
of the X-ray source, since the usual rapid variability outside the dip is absent during the 
obscuration. 

Of course, it is possible that the dip is due to intrinsic variation of the X-ray source. 
However, the origin of X-ray variability in AGN is not understood. Models which come close 
to successfully reproducing the observable quantities obtained from AGN lightcurves are of 
the shot-noise or 'rotating hot-spot' variety (e.g. Green et al. 1993; Bao and Abramowicz 
1996 and references therein). However, the parameters of such models must be highly tuned 
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in order to reproduce AGN power spectra. On the other hand, we show in this paper that a 
very simple-minded model can account for the temporal profile of the intensity dip in MCG 
—6—30—15 and briefly discuss the implications for AGN X-ray variability in general. 

2. A Simple Model for the Dip Temporal Profile 

Here we present a simple occultation model of the intensity dip (Figure 2). The dip is 
clearly divisible into several distinct time intervals. The turning points are labelled ti-tg in 
Figure 2 and the values shown for t 2 -t 7 are the nearest centers of the 512 s bins. The values 
for t\ and t$ are obtained from a model described below. Without any further analysis, 
one can make some immediate deductions, since the obscurer cannot travel faster than 
the speed of light. The duration of the observed dip (~3x 10 4 s) implies a source size 
< 9 x 10 14 cm and the duration of the absolute minimum, ~ 3000 s, implies that distance 
scales < 9 x 10 13 cm are resolved. 

The fact that there are two ingresses and egresses in the dip implies that either the 
obscurer or the X-ray emission must be spatially non-uniform. A non-uniform obscurer must 
have an extremely contrived shape and size relative to the X-ray source. For a non-uniform 
source, we can deduce a one-dimensional intensity profile. This profile must peak in three 
places to produce the nearly symmetric occultation profile observed. Two of the peaks must 
have similar intensity and be approximately equidistant from the central peak to preserve 
approximate symmetry. The obscurer must be larger than the distance spanned by the 
three peaks to avoid local maxima in the occultation profile. As the obscurer moves over 
the three peaks, if the first and last peaks are the more intense, then the size of the obscurer 
must be fine-tuned so that the last peak is covered at the same rate as the first peak is 
uncovered, in order to avoid a large, sharp rise after the absolute minimum. Moreover, the 
duration of the absolute minimum must be significantly greater than the duration of the 
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inflexions in the dip (£2^3 and t 6 -t 7 ), which is not observed. Therefore, the central peak 
must have the greatest intensity. Thus, the one-dimensional intensity profile must have the 
general (but not necessarily exact) form shown in Figure 2b. In this case, the first drop 
in intensity (t < £1) is small and not observable since most of the X-ray emission is still 
unobscured. The large drop over ti~t 2 is then due to the central source being covered. Note 
that the high data point when the source has just come out of the dip is not required to be 
explained by the model. This is because the central source, which we know exhibits rapid 
and erratic variability, is then completely uncovered. 

Although any two-dimensional model which gives the profile in Figure 2b is valid, it 
appears that an obvious and the least contrived realization of the inferred intensity profile 
is an X-ray source consisting of a high-luminosity central part, surrounded by a ring of 
dimmer emission. The gap between the central source and the ring is required, in order to 
fit the intensity profile between t 2 -t 3 and t 6 -t 7 . We construct a simple two-dimensional 
model which can explain the data (Figure 2b). The central source is represented by a circle 
of radius r s , and the obscurer by a circle of radius r$. In reality the latter may be more 
elongated in the direction of travel, but the important quantity is its length, 2r . The inner 
and outer ring radii are r\ and r 2 respectively. To account for the slight asymmetry in the 
dip (luminosity during t 6 -t 7 higher than t 2 -t 3 ), the half-ring uncovered last is made more 
luminous. Let the intensities per unit area of the central source, first and last covered halves 
of the disk be I s , h, I T respectively and let v be the transverse velocity of the obscurer. 
We find a fiducial model of the dip profile by adjusting r s /v, I s /h and h/h- The other 
parameters are fixed by the relations: (r x — r s )/v = £3 — t 2 , (r 2 — r-\)/v = £4 — £3 and 
2(r — r 2 )/v = £5 — £ 4 . A good fit (solid line in Figure 2a) is obtained with r s /v = 1536s, 
n/r a = 3.66, r 2 /r s = 7.32, r /r s = 8.32, I b /I s = .024 and I r /I s = .012. The dip slope 
between t\-t 2 strongly constrains the uncertainty in r s /v to < 20% of the fiducial value. 
The ranges on r±, r 2 and r depend on the uncertainties in r s and the turning points, and 
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again are < 20%, assuming a tolerance of 512 s (i.e. the bin width) for the turning points. 
The tolerances on h/h and I r /I s are also less than 20%. 

One interpretation is to identify the emission ring with the inner edge of an accretion 
disk. An inclined disk would explain the asymmetry in the occultation since Doppler effects 
due to rotation would give intensity enhancement or reduction from blueshifts or redshifts 
respectively. h/I r ~ 2 can easily be obtained for velocities of only 0.17c (Reynolds and 
Fabian 1997). In practice one must integrate Keplerian velocities over the disk (note, at 
6r g , v = 0.4 c). The persistent continuum emission during the dip is likely to be due to any 
uncovered part of the accretion ring plus emission from the rest of the accretion disk (which 
must have much lower surface emissivity than the inner ring). 

3. The X-ray Continuum and Fe-K Emission Line During the Dip 

We have shown that a simple-minded occultation model successfully reproduces the 
observed intensity dip in MCG —6—30—15. WY98 showed that the same model successfully 
explains the shape of the Fe-K line during the dip. The peculiar shape during the dip 
(bloated red wing and diminished blue-side emission) has been interpreted by a number 
of authors as evidence for line-emission from within 6r g of black hole (196; Reynolds and 
Begelman 1997; Dabrowski et al. 1997). In the occultation model, the shape is due to 
obscuration of the most blueshifted and redshifted part of the line emission from the 
putative accretion disk, leaving only redshifted emission. Line emission from within 6r g is 
not required. Moreover, WY98 find, in contrast to 196, that a huge Fe-K line equivalent 
width (EW > 1000 eV) during the dip is not required but is of the same order as the EW 
from the 4.2-day time-averaged value (~ 400 eV). Further, the Fe-K line intensity is not 
required to increase during the dip, compared to the time-averaged or flaring-state value (in 
contrast to the 196 analysis). In the analysis of WY98 the data are consistent with the line 
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intensity decreasing during the dip. 

WY98 also re-analysed the X-ray continuum spectrum during the dip. 196 interpreted 
the change in the 3-10 keV continuum power-law photon index, V, from ~ 2 (time-averaged 
'normal' value) to ~ 1.7 (dip value), as an intrinsic continuum change. WY98 found that 
the dip data are consistent with T fixed at the 'normal' value of 1.98, but allowing for 
extra absorption (equivalent neutral Hydrogen column density of ~ 2 x 10 22 cm -2 ). WY98 
interpreted this as the increase in opacity of the photo-ionized absorber in MCG —6—30—15 
(e.g. Reynolds et al. 1995; Otani et al. 1996) due to the central ionizing source being 
blocked by the putative optically-thick obscurer. Thus the occultation model for the dip 
temporal profile is entirely consistent with the observed Fe-K line and X-ray continuum 
during the dip. 



If accreting matter is exposed to the same UV/X-ray luminosity that we observe 
(L ~ 4 x 10 43 ergs s _1 ) then for gravitational infall to overcome outward radiation pressure 
requires the mass of the central black hole, M B h, to exceed 3 x 1O 5 M . Thus, if r\ is 
the radius of the inner edge of the accretion disk, identified with the last stable orbit of 
matter, and At is the time taken for the obscurer to traverse r 1 (i.e. t 3 — 0.5[t 2 + h]), 
then Kr g = ri where k = 6 or 1.24 for Schwarzschild or maximally rotating Kerr metrics 
respectively. But, n < cAt and r g = 1.48 x 10 13 (M BH /M ) cm, so M BH < 1.9 x 1O 8 M , 
or < 9.1 x 1O 8 M , for a Schwarzschild or Kerr black hole respectively. Assuming instead, 
the Keplerian velocity at the inner disk edge as the maximum velocity of the obscurer 



(v/c = \jr g /r = yl//c) yields smaller mass upper limits of M B h < 7.7 x 1O 7 M and 

-^bh < 8.2 x 1O 8 M for a Schwarzschild and extremal Kerr metric respectively. If 



the obscurer is at a distance, d, from the central source, r\ = vAt = cJr g /d gives 



4. Discussion 





- 9- 



d = c 2 [At] 2 /K 2 r g = 6 x lO 15 [At] 2 /[/t 2 (M BH /M )] cm. Thus, M BH > 3 x 1O 5 M implies 
d < 1.8 x 10 16 cm and d < 4.2 x 10 17 cm for Schwarzschild and Kerr metrics respectively. 
The origin of the optically-thick blobs is unspecified but Guilbert and Rees (1988) presented 
some simple arguments for the existence of dense (n > 10 15 cm -3 ), optically thick matter 
residing at the heart of accreting sources. The only independent estimate of the size of the 
'blobs' is that they should be much thicker than 10 9 /ni5 cm ( n i5 m units of 10 15 cm -3 ). 
The blobs must be optically thick even near their physical boundaries (i.e they must have 
fairly sharp edges), otherwise the dip profile would not be so well defined. Also, the blobs 
must be fairly stable, esepcially if they are created in the central region itself and 'propelled' 
up to high altitudes. 

The nature of the bright central X-ray source is intriguing. An inclined jet is unlikely, 
since even at 30° the inflexions in the dip profile (t 2 -t3 and t 6 -t 7 ) would have different 
durations. If the central X-ray source extracts its energy directly from the black hole then 
the metric is likely to be Kerr since energy cannot be extracted from a non-rotating black 
hole (Blandford and Znajek 1977; see also Ghosh and Abramowicz 1997). Occultations 
such as the one described here may occur frequently in AGN, but the relative sizes of the 
obscurer and source must be just right in order to observe such a clear event. Indeed the 
usual rapid variability or flicker, may be partly caused by the transit of optically thick 
bodies smaller than the source. We tested this hypothesis, again using a simple-minded 
model. Representing the bright central source as a circular disk with uniform emissivity 
(ignoring the ring due to its weaker emission), optically-thick blobs with ranges in radii 
(relative to the source) and velocities taken from Gaussian distributions were passed over 
the source. An additional parameter is required to specify the 'blob birth-rate' (i.e. the rate 
at which new blob trajectories are started). Such a model was used to produce predicted 
lightcurves for different model parameters. The power spectrum of each lightcurve was 
computed using the method of Papadakis and Lawerence (1993), omitting Poisson noise. 
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In the range ~ 10~ 5 Hz to ~ 10~ 2 Hz, no preferred or 'universal' power-law spectral slope 
was found. It is possible to produce power-law spectra with slopes similar to those typically 
measured (~ —1 to —2) but for most parameter values, the slopes are too steep. Thus, 
fine-tuning would be necessary to explain the 'universal' power-law slopes found in the 
handful of AGN in which it can be measured (e.g. Lawerence and Papadakis 1993; Green 
et al. 1993). This is essentially because if the product of blob crossing-time and birth-rate 
is too large or too small, there will be no variability. The direct simulated lightcurves 
assume, of course that the intrinsic source intensity is constant, which almost certainly is 
not the case. Thus our model does not explain AGN variability in general but its effects are 
potentially important to consider in any model of AGN variability. 

We thank the ASCA TEAM and mission operations at ISAS, Japan, for their efforts 
and hard work; Kim Weaver for her work on the dip spectrum, and Karen Leighly, Paul 
Nandra for some useful discussions. We also thank the anonymous referee. This research 
made use of archival data at the HEASARC, Laboratory for High Energy Astrophysics, 
NASA/Goddard Space Flight Center. 
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Figure Captions 



Figure 1 0.5-10 keV X-ray lightcurve from MCG-6-30-15 obtained by one of the CCD 
detectors (known as SISO) aboard ASCA (Tanaka et al. 1994). The data are binned at 512 
s. The reference time is 23 July 1994 UT 05:56:13. 

Figure 2 (a) Close-up of the latter portion of the lightcurve in Figure 1. Dashed line 
is the simultaneous background (not subtracted from the on-source data). The attitude 
of the satellite was stable to better then 0.004° and the apparent occultation event is 
consistent in all four ASCA detectors. Solid line is the predicted occultation profile of a 
simple two-dimensional model consisting of central source (radius r s ), and ring with inner 
and outer radii, r\ and r 2 respectively. The intensities per unit area of the central source, 
first-covered and last-covered halves of the disk are I s , h, and I r respectively. The obscuring 
body has a radius r . However it need not be circular, the important quantity is the length 
in the direction of travel. The solid line corresponds to a model with I s /h = 0.024 and 
h/ Ir = 0.012, r s /v = 1536 s, where v is the transverse velocity of the obscurer. The other 
parameters, r\jr s = 3.66, r 2 /r s = 7.32, and r /r s = 8.32, are determined by the turning 
points, £2^7 which were estimated from the data to the nearest 512 s bins. The time 
intervals shown correspond to the intervals between the front-end of the obscurer moving 
over the regions marked by the corresponding dotted lines (see text), (b) Schematic of the 
source geometry (top) and one-dimensional intensity profile (bottom) of the ring system, 
summed in a direction perpendicular to the velocity of the obscurer. In this model, the 
center of the obscurer passes over the center of the source but this need not be the case. 
Note that the high data point when the source has just come out of the dip is not required 
to be explained by the model. This is because the central source, which we know exhibits 
rapid and erratic variability, is then completely uncovered. 
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